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Q ■ Abstract For a skew normal random sequence, convergence rates of the distribution 

ly-^ I of its partial maximum to the Gumbel extreme value distribution are derived. The 

asymptotic expansion of the distribution of the normalized maximum is given under 
'■ an optimal choice of norming constants. We find that the optimal convergence rate of 

the normalized maximum to the Gumbel extreme value distribution is proportional to 
1/ log n. 
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^ : 1 Introduction 



fSJ ■ The biggest weakness of the normal distribution is its inability to model skewed data. This has 

led to several skewed extensions of the normal distribution. The most popular and the most studied 
of these extensions is the skew normal distribution due to Azzalini (1985). A random variable X 
is said to have a standard skew-normal distribution with shape parameter A G i? (written as 
X ~ SN(A)) if its probability density function (pdf) is 



(N 



X 



f\{x) = 2<j){x)^{Xx), -oo < X < +00, 

where (/»(•) denotes the standard normal pdf and <&(•) denotes the standard normal cumulative 
distribution function (cdf). It is known that SN(0) is a standard normal random variable. 

The skew normal distribution has received more applications than any other extension of the 
normal distribution. Its applications are too many to list. Some applications of the skew normal 
distribution that have appeared in the past year alone include: the distribution of threshold voltage 
degradation in nanoscale transistors by using reaction-diffusion and percolation theory (Islam and 
Alam, 2011); population structure of Schima superba in Qingliangfeng National Nature Reserve 
(Liu et al., 2011); rain height models to predict fading due to wet snow on terrestrial links (Paulson 
and Al-Mreri, 2011); modeling of seasonal rainfall in Africa (Siebert and Ward, 2011); modeling of 
HIV viral loads (Bandyopadhyay et al., 2012); multisite flooding hazard assessment in the Upper 
Mississippi River (Ghizzoni et al., 2012); modeling of diabetic macular Edema data (Mansourian 
et al., 2012); risks of macroeconomic forecasts (Pinheiro and Esteves, 2012); modeling of current 
account balance data (Saez et al., 2012); automated neonatal EEG classification (Temko et al., 
2012). 



1 



The aim of this note is to estabhsh the convergence rate of the distribution of the maxima for 
samples obeying SN(A). Chang and Genton (2007) showed that SN(A) belongs to the domain of 
attraction of the Gumbel extreme value distribution A{x) = exp(— exp(— x)). Rates of convergence 
of the distribution of maxima for a sequence of independent SN{0) random variables were studied 
by Hall (1979), Leadbetter et al. (1983) and Nair (1981). Precisely speaking, Leadbetter et al. 
(1983) proved that 

e-^exp(-e-^) (log log n)^ 

$ [anx + Pn) - A(x) ~ — 

16 log n 

holds for large n with normalized constants a„ and given by 

a„ = (21ogn)"^ and /?„ = a;;^ - ^ (loglogn + log47r) . 

The optimal uniform convergence rate of <I>"(a„x + 6„) to A(x) due to Hall (1979) is 

C2 



< sup 

log n ^giR 



(anx + bn) - A{x) 



< 



log n 



for some absolute constants < Ci < C2 with normalized constants a„ and 6„ determined by 

27r6^ exp ^5^^ = n^, = b~^. 
The following more informative result was established by Nair (1981): 

bllbl (^" {anX + K) - A{x)) - k{x)A{x)^ (u]{x) + ^^y^) A(rE), 

where the normalized constants a„ and 6„ are given by 

1 - $ (6n) = a„ = 6„\ 

where k{x) and iJ{x) are defined as 

■K{x) = 2'^ {x^ + 2x) e'"^ and uj{x) = -8"^ {x"^ + 4x^ + 8x^ + 16x) e'"" . 

The contents of this note are organized as follows. Section [2] derives some preliminary results 
related to SN(A) like Mills inequalities. Mills ratios, and the distributional tail representation of 
SN(A) for A 7^ 0. Convergence rates of the distribution of the maxima of SN(A) and related proofs 
are given in Section [3l In the sequel we shall assume that the shape parameter A / 0. 



2 Preliminary results 

In this section, some preliminary but important properties about SN(A) are derived. These 
properties not only imply that SN(A) belongs to the max-domain of attraction of the Gumbel 
extreme value distribution but they also help us to find two pairs of norming constants. 

The following Mills inequality and Mills ratio about SN(0) due to Mills (1926) are needed in 
this section, i.e., 

x"^ (1 + x-2)"^ (/)(x) < 1 - «>(x) < x^V(2;) (2.1) 
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for all X > and 

1 - <^{x) 1 



[X] X 



(2.2) 



as X —7- oo. For some improved Mills inequalities, see Mitrinovic and Vasic (1970, pages 177-180) 
and references therein. 

First we derive Mills inequalities and Mills ratios of SN(A), which are stated as follows. 

Proposition 1. Let F\{x) and fx{x) denote the cdf and the pdf of SN(A), respectively. For all 
X > 0, we have 



(i). ifX>0, 



which implies 



as X ^ oo; 
(a), if \ <0, 



which implies 



as X ^ oo. 



Proof. For x > 0, 



'1 /'i I „-2 



X (1 + x ) <-^^<- ' (2.3) 



1-Fa(x) 1 
f\{x) X 



< 



1 - F^^x) 
fx{x) 



l-Fx{x) 1 
/a(x) '^(1 + A2)x 



(2.4) 



(2.6) 



I f t^\ 

-2 I <^>(At)exp^--jdt 

1 f t^\ 

> J -^{Xt)expi^--jdt 

x"^^>(Ax) exp (^-y^ - J ^i^i) exp 
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Hence, 



> .-'<i.(A.) (4) + r *" (-H^'O 

In the case of A > 0, by (|2.8p . we can get 

(l + ^) ^ ^(AO exp (^-^^ > x-i$(Ax) exp (^-y^ , 

i.e., 

1^ $(At) exp (^-|^ > (^1 + $(ArE) exp (^-y) . 

By using (|2.7p and (j2.ip . we have 

^ $(At)exp (^-^^ 

< .-c|>(A.) exp (4) + f exp [-'-^t^) dt 
<x-icD(Ax)exp(-yj 1+^ ^ ' 



$(Ax) exp ( 



<.-..(A.)e.p(4)(,-*^)-\ 



Combining (fZO]) with (jZTOj) . we can derive (jZS]) and ([231). 
In the case of A < 0, by (|2.8j) and p.lh . we have 



l + $(At)exp(^-^ )dt 

„2 



>x-i<D(Ax)exp(-yj 1+ 



(l-«>(-Ax))exp(^-f 
>x-<^(Ax)exp(4) + 



For X > 0, 



[ HXt) exp (4) < x-c^(Ax) exp (4) (l " ^A^ (l + (t4)^) 



Combining ([2TT]) with ([232]) . we can derive ([23]) and ^^6^. □ 

By Proposition [H we can derive the distributional tail representation of the skew normal distri- 
bution. This representation is useful to find optimal normalized constants to establish expansions 
of the distribution of the maxima for SN(A) samples. Similar expansions for SN(0) samples were 
given in Nair (1981). 

Proposition 2. Let F\{x) and f\{x) denote the cdf and the pdf of SN {X) , respectively. Then, 

l-FA(x) = c(x)exp (^-j^^^^dt^ 

for large x, where c{x), g{x) and the auxiliary function f{x) are determined according to the sign 
of X, i.e., 

(i). for A > 0, 




as rr — )• oo, 



f(x) = - > with fix) 
x 



— )■ as X — )• cx) 



(2.13) 



and 



g{x) = 1H — ^ ^ ^ OLS X ^ oo. 



(2.14) 



(a), for A < 0, 



expi-i±^ 



fix) = ,^ /.^s > with fix) = - \ ^ as x ^ oo (2.15) 
(1 + A^)x (l + A^jx"' 

and 

2 

9ix) = ^ + (i + ;^2)^2 ^ ^ x^oo. (2.16) 

Proposition [2] shows that Fx G D{A) by Corollary 1.7 of Resnick (1987). The norming constants 
a„ and 6„ can be determined by 

l-Fxibn) = l/n, a„ = /(6„) (2.17) 

such that 

lim Fj:(a„x + 6„) = A(x). (2.18) 

n— >oo 

By using Mills ratio of the skew normal distribution and arguments similar to the case of SN(0) 
provided in Leadbetter et al. (1983), we can choose another pair of normalized constants such that 
(f2T8]l holds. 
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Proposition 3. Let {Xn,n > 1) be a sequence of independent and identically distributed random 
variables with common cdf Fx{x). Let = max(Xfc,l < k < n) denote the partial maximum. 
Then 

lim P (M„ < anx + Pn) = A(x), (2.19) 
where the are given by 



2(2 log n) 2 



for A > 0, and by 



^^ , \2\-5/oi R /21ogn\2 loglog?i + log(-27rA 

Vl + ^V (1 + A2)2 (21ogn)2 



for A < 0. 

Proof. First we consider the case of A > 0. As -Fx is continuous, for x £ R, tliere exists Un = Un{x) 
such that n(l — Fx{un)) = e~^. By using ()2.4p . we have 

as n — 7- oo, and so 

1 

log n - log Un + - (log 2 - log vr) - ^ + log ^> (Atx„) + x ^ (2.20) 
as n — )• oo, which implies 



2 log n 

as n — )• oo. Taking logarithms, we have 

2 log Un — log 2 — log log n — )• 

as n — )• oo, so 

logun = i(log2 + loglogn) + o(l). (2.21) 

Noting that hm„^oo ^{^Un) = 1 for A > and putting (|2.2ip into (|2.20p . we have 

/ 1 loglogn + logvr x f 1 W ^ 

Un = (2 log n) 2 1 — h h O 



2 log n log n \ log n 

a A loglogn + logvr x / 1 

(21ogn)2 1 — — + -^^ + o' 



4 log n 2 log n \ log n 

^ , ^oi \h loglogn + logvr / 1 

Y + (2 log n) 2 ^ — + o r 

(2 log n) 2 2(2 log n) 2 y(logn)2 
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= UnX + /3n + (an) ■ 

So, by Theorem 1.2.3 and Theorem 1.5.1 of Leadbetter et al. (1983), we have 

hm P (M„ < anx + /3„) = A(x) 

n— >oo 

with 

an = (21og?i)~^ 

and 

« ^n^ \h loglogn + logvr 
/3„ = (2 log n) 2 - 



2(2 log n) 



Similarly, for A < there exists Vn = Vnix) such that n{l — Fx{vn)) = e ^ . By using \2M and 
2.2p . we have 



nv-' (1 + A^)-^ (-Avr)-i exp [x - ^u^) ^ 1 

as n — 7- oo. By arguments similar to the case of A > 0, we can derive 

lim P (M„ < anX + /3„) = A(x) 

with 

an = (l + A^)"^ (21ogn)-^ 

and 

^ /21ogn\2 loglogn + log(— 2A7r) 
Vl + AV (1 + A2)5 (21ogn)5 

The proof is complete. □ 

3 Convergence rates of extremes 

In this section, we establish two different convergence rates of the distribution of M„: one for 
the norming constants a„ and /3„ given by Proposition O and the other for the norming constants 
a„ and 6„ determined by (|2.17p . 

Theorem 1. For norming constants an and /3n given by Proposition {31 we have 
(i). For A > 0, 

/ n \ A / N A(2;)e~^ (loglogn)2 

F^ {anX + I3n) - A{x) ^ \; V s s ; 

lb log n 

as n ^ oo; 
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(ii). For A < 0, 



^r, , ^ N . / N A(x)e ^ (log log n)^ 

4 log n 



as n oo. 



Proof. Set tin = anX + /3n and = n(l — F\{un)), where q„ and are given by Proposition [3l 
For A > 0, 



implies 



n.„ = (21ogn) 2x + (21ogn)2 — 2 -'^(21ogn) 2 (log log n + log vr) 



2 , log log n + log vr (loglogn)2 

u„ = 2 log n + X + (1 + 0(1 ) 

' z lb log n 



,-,2.o.„,4(:.o(»)) 



and 



Noting that log<^(Ati„) = o ((log log n)^/ log n) for large n and rewriting ()2.3p as 

i fl _ < l--^A(:g) < i + 'A(Ax) 

x\ l + x^y /A(a;) a; \ Ax — (^(Ax) 

we have 

= 77, (1 - Fa (n„)) 

2n / u. 



exp (-^ + log <^ (An„)) (l + O (u^^)) 

(i+'^a))^ 



27rn.^ 

^ _ (log log 77)^ 
16 log 77 



Obviously, for t(x) = e 



/ N / N -x(loglogn)2 „^ (log log 77)2 
t[x) — Tn{x) = e (1 + o(l)) ~ e ^ 

16 log 77 16 log 77 

for large 77. So, by using Theorem 2.4.2 of Leadbetter et al. (1983), we can get (j3.ip . 
In the case of A < 0, using (|2.5p and the following fact 



i-Hx) = ^[i + o{x-^)], 

we can get (j3.2p by arguments similar to the case of A > 0. □ 
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Theorem [2] shows that the convergence rates of the distribution of the normahzed maxima are 
different even though both (j2.18p and (|2.19p hold, which imphes a„/a„ — )■ 1 and (/3„ — &n)/on — ^ 
by Khintchine Theorem (cf. Leadbetter et al. (1983), Resnick (1987)). The auxihary function f{x) 
of the distributional tail representation of Fx(x) plays an important role on the convergence rates 
for distributions belonging to the domain of attraction of the Gumbel extreme value distribution. 
For related work, see Peng et al. (2010) and Liao and Peng (2012): Peng et al. (2010) derived 
convergence rates of the distribution of maxima for random sequences obeying the general error 
distribution; Liao and Peng (2012) derived convergence rates of the distribution of maxima for ran- 
dom sequences obeying the lognormal distribution. Another interesting application of the auxiliary 
function f(x) is to characterize F G D{A). For more details, see Goldie and Resnick (1988). 



Theorem 2. For norming constants a„ and given by (|2.17p . we have 



bl Fl {anx + bn) - A(x) - k(x)A(x) 



oj{x) + 



K?{x) 



(3.3) 



as n ^ oo, where both k{x) and uj{x) may depend on A according to the sign of X as: 
(i). for A > 0, 

k(x) = 2"^ (x^ + 2x) e-"", uj{x) = -8"^ {x'^ + 4x^ + Sx"^ + 16a;) e""^; 
(a), for A < 0, 

k{x) = 2-1 (1 + A^)"^ (x^ + 4x) e"^ 

and 



uj(x) 



l + X^) ^ \x^x* + 8A^x^ + 24A^x=^ + 16 (l + 3A^) 



To prove Theorem [2l we need two lemmas. The first lemma is about a decomposition of the 
distributional tail representation of the skew normal distribution. 

Lemma 1. Let F\{x) denote the cdf of the skew normal distribution. For large x, we have 
(i). for A > 0, 



1 - Fx{x) 



Ivre / 



1 - + 3x-'^ + O (x"^) $(Ax)exp 



fit) 



dt 



(3.4) 



with f{t) and g{t) given by (|2.13p and (j2.14p . respectively; 
(a), for A < 0, 



l-Fx{x) 



exp 



1 



1 + 3A2 



-X 



(_A)vr(l + A2)L A2(l + A2)- 

x~^ + 0(x" 



15A'^ + 10A2 + 3 _4 



A4(l + A2)2 

with f{t) and g{t) given by (j2.15p and (j2.16p . respectively. 



exp 



fit)' 



dt (3.5) 
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Proof. By integration by parts, we have 



1 - Fxix) 



1 - + Sx""^ - 15x-^ + 



A H^x) _ A (3 + A^) ^{Xx) 
(1 + A2)2 $(Ax) 



A 



3 + 



1 + A2 ^>(Ax) 
4(3 + A2)\ ,/,(A, 



1 + A2 I ^ (1 + A2)2 i «^(Ax) 



15 + 



3 + 



4(3 + A2)\\ (t){\x) _j 



1 + A21 ^ 1 + A2 I ^ (i + A2)2 M $(Ax) 



210 o , / t2\ 



+ 



8A 



7r(l + A2) 



18 \ 



+ 



192A(3 + A2) 
^(1 + A2f 

It is easy to check that both 



t exp 



dt 

t exp 
1 + A2 



1 + A2 



dt 



e dt. 



(3.6) 



< / t^^$(At)exp 
and 



t' 



dt < X ^$(Ax) exp + 



27r Jx 



■00 



r exp 



1 + A 



2 i'jdi(3.7) 



< / t"^ exp - 



1 + A^ 
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1 + A2 



exp 



1 + A2 



(3.8) 



hold for X > 0. 

For A > 0, it is obvious that x^4'{Xx) — )• and ^{Xx) — )• 1 as 2; — )• 00. So, by (|3.6p . (|3.7p and 
(j3.8|) . we can get 

= ' (^1 - + 3x-^ + O {x-^) ) $(Ax) exp 

which is (jSli]) . 

In the case of A < 0, first notice that 

ct){x) 



t(l + ^]dt 



1 _ <I>(^) = ^ [1 _ ^-2 + 3^-4 + o (^-6)] 



(3.9) 



for large x, cf. Castro (1987). So, by (|3.6p and p.9p . we have 



1 - Fx(x] 



2(l){x) 



$(Ax) (1 - + 3x-4 + O (x-6)) + 



A (/)(Ax) A (3 + A2) 0(Ax) 



(1 + A2)2 x= 



+ 



A 

TTa2 



3 + 



4(3 + A2)\ (j){Xx) 



(1 + 



A2)^J 



+ 



0(Ax) 
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exp^-i^x^J . 1 + 3A^ 15A^ + 10A^ + 3 

(_A)vr(l + A2)x2 V A2(l + A2)'^ + A^ (1 + A^)^ ^ ^ ^ 
exp(-ii^) ^^3^2 _^ 15A4 + 10A2+3 .4 



1 — --T7- r^X H r5 X 



(_A)vr(l + A2)V^ A2(l + A2)-^ ' (1 + A2)2 

+0 {x-') ) exp (- (1 + A^) t (1 + (Y^^) rfi) , 

which is ()3.5p . The proof is complete. □ 
To prove (j3.3p . we need the following auxiliary result. 

Lemma 2. Lei H\{hn',x) = F\{anX + bn) and hx{bn',x) = nlogHx{bn;x) + e~^ with the normalized 
constants an and bn given by (j2.17p . Then, 

lim (6„;x) - ) =uj{x), (3.10) 

where k{x) and uj[x) are those given by Theorem\^ 

Proof. First, we consider the case of A > 0. Obviously, 6,„ — t- 00 if and only if n — )• 00 since 
1 — Fxibn) = n^^. The following two facts hold by Proposition [TJ 

Um <'-^^<^-^ (3.11) 

/A {bn) 



and 



1 - Fx [bn + b~^x] 
lim ^^^3 = 0. (3.12) 

bn 



For simplicity, let 

a>(AM (1-6-2 + 36-4 + 0(6-6)) 



^{\{bn+bn^x))(l-{bn + bn^x) ^ + 2, {bn + bn^ x) Vo((fe„ + fen^x) 



Then limn-s>oo ^A(&n) = 1 and 
^A {hn) - 1 = (1 + 0(1)) {Xbn) - (A (6„ + b-'x)) + 6^^ {q> (A (6„ + 6;^x)) - <I> (A6„)) 
-2xb-^<^> (A {bn + b-'x)) + 36;;^ (A6„) - $ (A (6„ + b~'x)) (l + 6;:2^)"' 
+0 (6-6 



So, we can check that both 



and 



lim = (3.13) 



hm = _2x (3.14) 



11 



hold. By (|3.4p . we have 

1 - Fx (bn) 



1 - Fa {bn + bn^x) 



Axibn)exp / {bjy + 



1 



A 



1 

2 



1+ / y + 



K y + 



bl (1 + bn^y) 
bl{l + bn^y) 



]dy^ 
dy 



(l + 6n'y) 

Combining ([23]), (fSlTT) . (|3T^ . (|3T^ . (IXTil) and (IXTSD together, we have 

logi^A (6n;a;) + n~^e"^ 



(1 + 0(1)) . (3.15) 



hm bnhx{bn;x) 



Hm 

n— >oo 



n 



"71 



Hm 

n— ^-oo 



log Fa (6n + b-\x)) + [1 - Fa (6^)] e"" 

/a (bn) bn"" 

- [1-Fa(6„ + 6;:1x)] -i [1-Fa(6„ + 6^1x)]'(1 + o(1)) 



lim 

n— s-oo 

f\{bn)bn 



fx{bn)bn 



l-n(bn) 



lim 

n— ^-oo 



bn [l - Fa (6n + 6„ ^a;)] l-FA(6„+b;r'a:) 



fx{bn 



bn' 



e ^ hm 



Ax {bn) -l + Ax {bn) /o" {bnh + ,.(,^\-.^) ) (1 + o(l)) 



e lim 

^1 1^2 



^^l + bn'y 
(.x2 + 2x) e""^ = k{x), 



dy 



(3.16) 



where the last step follows by the dominated convergence theorem. It remains to prove (j3.10p for 
A > 0. By arguments similar to those of (|3.16p . we have 



lim b^ 

n^oo 



bljix {bn;x) - k{x) 



lim 

n— >oo 



lim 

n— >oo 



2x^+2x 



lim 

n— >oo 



e ^ hm 



logHx{bn;x) + n e ^ — b^n k{x) 
log Fa {bn + b-^x) + [1 - Fa (6„)] (l - 6" 

/a (fen) 

bn[l-Fx{bn + b-^x)]l-F^{b,^+b-'x)'' V 2 J 



/a (fen) 
^A (fen) fe^, 

. , x2 + 2x r 

-^A (fen) - 



^^1 + fen^y 



dy 



x"^ + 2x 



2 7 



1 



1 + fen ?/. 
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, Ax - 1 



1 + 2x 



2 462 



2x e~ 



2 2 V 2 

-8"^ (x'^ + 4x^ + Sx^ + 16x) e""^ = a;(x). 



Second, we consider the case of A < 0. By (j2.6p . ()2.15p . and arguments similar to the case of 
A > 0, we can see that 



bn 



lim 

n— >oo 



l-FA(&n+(l+A2) h~^X 



and 



lim 

n— >oo 



hold. Setting 



Bx {bn) 



l-Fx(hn + {l + \Y^b~^: 
1 1+3A^ ^-2 I 15A-*+10A^+3 j,-4 , 



(3.17) 



(3.18) 



1 



1+3A2 



bn + {l + \^)'^bn^X 



A^(l+A^) 

we have lim„_j.oo Bx{bn) = 1 and 



I ISA^+lOA^+S 
^ A'i(l+A2)2 



6„, + (l + A2)-i6;^ix +0{bn'^) 



-4 



Sa (6n)-l = (1 + 0(1)) 



2 (l +3A2) ; 

" A2(l + A2)' 



which implies 



and 



lim 



i?A (6n) - 1 _ 2(1 + 3A2) 



A2(l + A2) 



2 • 



(3.19) 



(3.20) 



By p.Sp . we have 



1 - i^A (6n) 



l-i^A(6n + (l + A2)-l6-lx] 



(^n) exp 



(1 + A2)fe2 1 + A2^2 



l + (l + A2)-l6-2y 
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1 + 



y 



+ 



1 



(1 + A2)62 1 + A2,2 

2 

+ 



l + (l + A2)-^&-2y 
1 



(1 + A2)62 1 + A2 2 



1 + (1 + A2)-i hn^y 



dy\ (l + o(l)) 



. (3.21) 



Combining (g^SD, (|XT7D . (fXTHD . (|XTm . and ([5:^ together, we have 



lim b^hxibnjx) = lim 



l-Fx bn+ l + X^) 'K^X 



(1 + A2)-Va(M 



e ^ lim 



Bx {hn) -l + Bx {hn) Jo ( 



2/(1+A2) 



(l + A^) ^ (x2 + 4x)e- 



and 



hm 6^ 

n— s-oo 



hm ■ 

n— >oo 



l-n(fe„ + (l+A2)-ife-ix)^ 



1 _ /,-2_xf+42_ 
^ 2(1+A2) 



e hm 



(l + A2)-VA(6n) 



y 2 

+ 



+ 1 + A2 1 + (1 + A2)-i 6;; 



dy 



-Bx {hn) bl- 



'2(1 + A2) 



-5a {hn) 



+ 4x 



y 2 

+ 



1 



2(l + A2)7o V1 + A2 ' l + A2i + (i + A2)-i6^-2^, '"^^ 



+ ^5a (^^n) 



y 2 

+ 



1 



2(1 + 3A2) 



V1 + A2 l + A2i + (i + A2)-i6,;r2y 

1 



dy] (1 + 0(1)) 



y 2 



dy 



^ a;2 _^ 43. 



+ A2 i + A2i + (i + A2)-i6-2yy 2(1 + A2) 
1 / 2;2 + 4x ^ ^ 



(l+o(l)) 



A2(l + A2)2 (1 + A2)2 2V2(1 + A2)^ 

= -8-1 A-2 (1 + A^)"^ [A^x^ + SA^x^ + 24A\2 + 16 (l + SA^) x] e"^ = w(x). 

The claimed result follows for A < 0. The proof is complete. 

Proof of Theorem [21 Obviously, by Lemma El we have hx{bn',x) — )■ and 



□ 



f^ fej-' {hn-.x) 

i=3 



< exp [hx {hn]x)] 1 



as n — >• 00. By using Lemma [2] once again, we have 



hi {F^ {anx + hn) - Hx)) - k{x)K{x) 



14 



bl (exp {hx x)) - 1) - k{x) A{x) 



ujix) + 



The desired result follows. 



bl {blhx (bn, x) - k{x)) + bf^hl (6„; x) [i^ + hx {bn, x) ^ 
k'^{x) 



°° hi-Hbn;x) 



i=3 



A{x) 



A(x). 



□ 



Theorem [2] establishes the asymptotic expansion of F^{anX + bn)- Meanwhile the convergence 
rate of F^{anX + 6„) to its limit distribution A{x) is proportional to 1/logn by Theorem [2] since 
one can check that l/bl = (1/logn) through (ETT]). 
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